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SUPPORTED  BY  GRANT  AFOSK-78-3662 
(September  30,  1979-Septomber  29,  1980) 

Basic  Instability  Mechanisms 
in  Chemically  Reacting  Subsonic  and  Supersonic  Flows 

T.  Y.  Toong  and  G.  E.  Abouseif 
Massachusetts  Institute  of  Technology 

Summary  of  Progress 

The  stability  of  one-dimensional  piston-supported  detonations  has 
been  examined  both  numerically  and  analytically.  Numerical  calculations 
were  conducted  for  a  one-step,  first-order,  irreversible  reaction  obeying 
an  Arrhenius  rate  expression.  An  approximate  linearized  stability  theory 
was  also  developed  for  the  case  of  high  act ivat ion-enerev  reactions  and 
the  mechanism  of  instability  identified. 

The  analysis  demonstrates  that  the  interaction  between  the  irreversibl 
temperature  fluctuations  and  the  reaction  zone  induces  an  oscillatory 
energy-source  field,  which  then  leads  to  shock  perturbations  and  thereby 
the  temperature  fluctuations.  The  frequencies  of  all  the  unstable  and 
stable  modes  of  this  system  can  be  predicted.  They  agree  remarkable  well 
with  the  findings  of  the  numerical  calculations  and  the  observations  it. 
blunt-bodv  flow  experiments. 

The  problem  of  direct  initiation  of  gaseous  detonations  was  also 
examined  theoretically,  to  determine  the  correlation  between  the  critical 


energy  and  power  for  successful  initiation,  and  to  predict  their  respective 
threshold  values.  The  results  were  found  to  agree  well  with  the  experi¬ 
mental  observations  of  cylindrical  detonations  in  oxy-acetylene  mixtures. 


One  paper  is  in  press  and  another  two  papers  have  been  prepared  for 
journal  publication. 
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I.  Objective  and  Scope  of  the  Work 


The  ultimate  objectives  of  this  work  are  three-roid.  First,  to 
determine  the  mechanism  of  development  of  the  transverse  wave  structure 
in  gaseous  detonations.  Second,  to  establish  a  means  of  es  t  i  ;nat  ing 
the  cell  size  of  self-sustained  waves.  Finally,  to  determine  the 
necessary  requirements  for  successful  direct  initiation  and  the  conditio 
leading  to  failure.  This  information  is  equally  important  for  the 
understanding  of  detonations  in  liquid  and  solid  explosives. 

In  the  first  part  of  this  work,  we  examine  the  phenomenon  of 
longitudinal  oscillations  observed  in  blunl-Dody  tlov:  experiments  and 
in  numerical  simulations  of  unsteady  one-dimensional  detonations.  Our 
main  objective  is  to  determine  the  physical  mechanism,  of  instability. 

In  future  work,  we  shall  extend  the  study  to  the  two-dimensional  case, 
to  analyze  the  transverse  wave  structure  of  detonation  waves. 

In  the  second  part  of  this  work,  we  examine  the  problem  of  direct 
initiation  of  gaseous  detonations,  to  determine  the  correlation  between 
the  igniter  critical  energy  and  average  power  for  successiul  initiation 
and  to  predict  their  respective  threshold  values. 

The  following  discussion  summarizes  the  results  of  these  two 


studies. 


II.  Results  and  Discussion 

(A)  Instability  in  Gaseous  Detonations 

Our  analysis  has  demonstrated  that  the  mechanism  of  instability  in 
gaseous  detonations  is  a  non-acoustic  one.  The  sequence  of  events  works 
out  as  follows:  The  pressure  waves  perturb  the  shock  front  and  generate 
irreversible  temperature  fluctuations.  The  latter,  in  turn,  are  carried 
by,  and  interact  with  the  reacting  flow,  thereby  inducing  a  fluctuating 
energy-source  field.  This  energy-source  field  then  generates  pressure 
waves  which  lead  to  shock  perturbations  and  thereby  the  temperature 
fluctuations . 

Figure  1  shows  the  temporal  development  of  the  shock  pressure  for  a 
one-dimensional  detonation  wave  at  Chapman-Jouguet  conditions.  One 
observes  several  unstable  modes.  The  first  (and  lowest-f requency )  mode 
has  a  period  of  %  20  characteristic  chemical  times,  and  a  peak-to-peak 
amplitude  of  %  4  times  the  steady-state  shock  pressure  p^  .  The  second 
mode  lias  a  period  of  x  1.5  chemical  times  and  a  peak-tc-peak  amplitude  of 
T.  0.25  the  p  .  The  figure  also  shows  higher  modes  of  shorter  periods 
and  smaller  amplitudes.  Preliminary  analysis  of  the  two-dimensional 
case  suggests  that  it  is  the  first  mode  which  produces  the  transverse 
wave  structure  in  mu  1 1 i -d imens iona 1  detonations. 

Figure  2  sh  rws  the  predictions  of  the  oscillation  modes  of  a 
typical  detonable  system  for  several  degrees  of  overdrive,  f.  The 
untagged  symbols  denote  the  results  of  the  present  theoretical  work  and 
the  tagged  ones  denote  the  results  of  an  exact  numerical  stability 
analysis  reported  in  the  literature.  u>.  and  are  the  imaginarv  and 
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real  frequencies  of  the  oscillatory  nodes,  ui.  >  0  denotes  an  unstable 
system  where  the  oscillation  grows  in  time.  The  figure  shows  that  the 
lowest-f requency  modes  are  the  most  unstable.  Moreover,  the  present 
theory  is  indeed  successful  in  predicting  the  periods  and  amplification 
(or  attenuation)  rates  of  all  the  oscillatory  modes. 

Figure  3  compares  the  predictions  of  oscillation  frequencies 
according  to  the  present  theoretical  work  (solid  lines)  with  the 
experimental  results  for  blunt-body  flows  in  hydrogen-oxygen  systems. 
Again,  the  present  theory  is  capable  of  predicting  the  periods  for 
both  the  high-  and  low-frequency  unstable  modes. 

Work  is  currently  in  progress  to  dot-ermine  the  unstable  frequencies 
for  two-dimensional  flows  and  to  predict  the  cell-siee  in  self-sustained 
waves . 

(B)  Direct  Initiation  of  Gaseous  Detonation 

A  simple  theoretical  model  was  developed  to  determine  the  correlation 
between  the  igniter's  critical  energy  and  power  fur  direct  initiation  of 
detonations  and  to  predict  their  respective  threshold  values.  The 
model  comprises  a  constant-velocity  piston  which  is  set  in  motion  at 
time  t  --  0 .  The  resulting  constant-velocity  shock  wave  heats  up  the 
explosive  mixture  and  triggers  the  chemical  reaction.  It  is  shown  that 
for  direct  initiation,  the  piston  should  remain  in  motion  ror  a  t i me 
interval,  at  least,  equal  to  the  induction  time.  The  predictions  based 
on  this  simple  model  are  in  good  qualitative  and  quantitative  agreement 
with  the  experimental  findings,  reported  in  the  literature  on  evlindrieni 
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detonations . 

Figure  4  shows  the  calculated  shock  pressure-time  history  for 

✓N  ^  ^  ~ 

different  intervals  of  piston  motion  t  .  For  t  <  t,  (where  t.  is 

P  Pi  i 

the  induction  time)  the  figure  shows  that  the  shock  decays  monotonically 

to  sonic  conditions.  For  t  >  t^,  the  shock  pressure  increases  very 

rapidly  and  then  decays  to  C-J  conditions.  For  t  =  t.,  the  figure 

p  l 

shows  an  initial  decay  in  shock  pressure  (similar  to  the  case  t  <  t.), 

P  i 

followed  by  two  consecutive  rapid  pressure  jumps.  Finally,  the  shock 
pressure  decays  asymptotically  to  the  C-J  value,  Pc_j  •  The  figure 
demonstrates  that  for  successful  initiation,  the  piston  should  remain 
in  motion  for  an  interval,  at  least,  equal  to  the  induction  time. 

Figure  5  shows  a  plot  of  critical  energy  versus  average  igniter 
power  P  for  an  oxy-acety lene  mixture.  The  solid  curve  denotes  the 
results  predicted  by  present  model,  while  the  symbol  (0)  denotes  experi¬ 
mental  results.  The  figure  demonstrates  that  the  present  model  is  indeed 
successful  in  predicting  both  the  qualitative  and  quantitative  behavior 
of  the  energy-power  requirements  for  direct  initiation. 

III.  Publications  and  Reports 

See  attached  Enclosure  1. 

IV.  Professional  Personnel 

Professors  T.  Y.  Toong  and  C.  E.  Abouseif 


V. 


Interact  ions 


Ref.  12  listed  in  Enclosure  I  was  presented  at  the  7th  International 
Colloquium  on  Gasdynamics  of  Explosions  and  Reactive  Systems,  University 
of  Gottingen,  Gottingen,  Germany,  August,  1979. 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig .  5 


CAPTION  TO  FIGURES 


Temporal  development  of  the  shock  pressure  for  a  one-dimensional 
deconation  wave  at  Chapman-Jouguet  conditions. 

Predictions  of  oscillation  modes  of  a  typical  detonable  system 
fur  several  degrees  of  overdrive.  Untagged  symbols  denote  the 
present  theory.  Tagged  symbols  denote  exact  solutions. 

Predicted  oscillation  frequencies  versus  observed  frequencies 
for  hydrogen-oxygen-diluent  mixtures.  The  solid  lines  indicate 
the  present  theory.  Symbols  denote  experimental  results. 

Calculated  shock  pressure-time  history  for  different  intervals 
of  piston  motion. 

Critical  energy  versus  average  igniter  power  for  an  oxy-acetylene 
mixture.  The  solid  curves  denotes  the  present  model, (0)  denotes 


experimental  results. 
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